Inelastic X-Ray Scattering Study of Exciton Properties in an Organic Molecular 

Crystal 
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Excitons in a complex organic molecular crystal were studied by inelastic x-ray scattering (IXS) 
for the first time. The dynamic dielectric response function is measured over a large momentum 
transfer region, from which an exciton dispersion of 130meV is observed. Semi-empirical quantum 
chemical calculations reproduce well the momentum dependence of the measured dynamic dielectric 
responses, and thus unambiguously indicate that the lowest Frenkel exciton is confined within a 
fraction of the complex molecule. Our results demonstrate that IXS is a powerful tool for studying 
excitons in complex organic molecular systems. Besides the energy position, the IXS spectra provide 
a stringent test on the validity of the theoretically calculated exciton wavefunctions. 

PACS numbers: 71.35.-y, 78.70.Ck, 31.15. Ct 
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In the past several decades, optical applications of or- 
ganic materials have been extensively explored. [l|, [U, [H, 
El IE IE Organic light emitting devices, photovoltaic 
cells, photochromic materials, biosensors, and nonlinear 
optical devices among many others have generated a lot 
of interests. Great effort has been invested on designing 
new functional materials with special optical properties. 

As the optical processes are largely dominated by exci- 
tonic excitations in these materials, the detailed informa- 
tion of excitons is crucial. Particularly, it would greatly 
benefit the material design if one knows where the exci- 
ton is located in a molecule, what its spatial extent is, 
and how the exciton energy and dispersion are associated 
with the local structure. 

These basic yet crucial questions on exciton cannot 
be addressed by conventional techniques such as Raman 
scattering and absorption, as they can only give infor- 
mation near zero momentum transfer. So far, high en- 
ergy electron-energy-loss spectroscopy (EELS) has been 
applied successfully to thin film samples of conjugated 
oligomer or poly mer such as a-6T, 6P, carotenoid, and 
TPD. @, i, [Hi [lj. EELS is limited to the low mo- 
mentum transfer q due to strong multiple scattering at 
high q, which would severely smear out the spectrum, 
and the diminishing matrix element (~ q~ 4 )- As a re- 
sult, EELS studies were mostly focused on planar n 
conjugated molecules where the excitons are more ex- 
tended in space, and information at low q is sufficient 
for understanding the exciton behaviors. On the other 
hand, the majority of organic optical materials are com- 
posed of small molecules with complex local structures. 



'Electronic address: dlfeng@fudan.edu.cn; zgshuai@iccas.ac.cn 



Even for polymers, the optically functional portions are 
mostly small sectors attached to the backbones for bet- 
ter mechanical performance. For small molecules or small 
molecular clusters, the size of the excitons is presumably 
very small, and high-q 1 information is crucial for under- 
standing their behaviors. The properties of excitons in 
these systems remain largely unexplored experimentally. 
Theoretically, although many quantum chemical meth- 
ods have been employed to study these systems, their 
feasibility lacks strong experimental support. 

Inelastic X-ray scattering (IXS) has proven to be a 
powerful tool for investigating the electronic excitations 
in inorganic systems. [12| For example, IXS data have 
helped our understanding of the metal-insulator transi- 
tion [fI], plasmon excitations [3], and band gap [l5j |. 
IXS is a clean and direct probe of the dynamic structure 
factor S(q,u>), which is proportional to q 2 x 7m(e _1 ), 
e being the dynamic dielectric function. In comparison 
with EELS, IXS is almost free from multiple scattering 
effects, and it can generate reliable information at high q. 
Therefore, IXS is particularly advantageous for studying 
excitons that are more constrained in space. 

In this paper, we report the first IXS measurement of 
the Py-SO organic molecular crystal. The dynamic di- 
electric responses of such complex small molecules were 
measured over a large momentum transfer region. The 
dispersion of the lowest exciton was observed to be 
about 130 meV. Quantum chemical calculations based 
on ZINDO (Zerner's intermediate neglect of differen- 
tial overlap) /SCI (single configuration interaction) [l6| 
was performed on both molecular aggregate and a sin- 
gle molecule. We found that the ZINDO/SCI method 
captured very well the internal structure of the excitons. 
As the momentum dependence of the dielectric function 
gives strong constraints on the exciton wavefunctions, the 
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FIG. 1: A unit cell of Py-SO molecular crystal. 



good agreement between the calculated IXS spectra and 
the experiments indicates unambiguously that the low- 
est Frenkel exciton distributes in a fraction of a single 
molecule. Our results provide comprehensive informa- 
tion for understanding the properties of excitons in such 
complex molecular systems. 

Spirooxazines (SO) are a type of well-known pho- 
tochromic compounds. Open-ring photomerocyanine 
form of spirooxazines (Py-SO, C21N4O2H22) was syn- 
thesized recently pjj, which exhibited alkali induced 
chromism and the thermochromism in alkali medium. 
High quality single crystal Py-SO samples, typically in 
the size of 0.8mm x 2mm x 150/im, were prepared by re- 
crystallization. Fig. [T] shows one unit cell, which contains 
two molecules. The bottom right Py-SO molecule, from 
left to right, consists of one phenyl ring, one aromatic 
ring, and a C-N bond bridge to the quinoid on the right. 
At room temperature, the triclinic lattice parameters are 
a = 8.356A, b = 9.510A, c = 12.096A, a = 89.93°, 
13 = 75.58°, 7 = 81.92°. 

IXS measurements were carried out in transmission 
mode at the Taiwan Beamline BL12XU at SPring-8.(l8j 
A Si(444) spherical analyzer with 2m radius of curvature 
was used. A Si(400) high-resolution monochromator was 
used to scan incident photon energy around 7908. 75eV. 
The total energy resolution was estimated to be about 
170meV based on the FWHM of the quasi-elastic lines of 
the sample. Momentum resolution was about 0.14A -1 . 
Data were taken at room temperature. Great care was 
taken to avoid beam damage to sample during the exper- 
iments by changing the probing spot every scan, since the 
beam spot was only about 120/im x 80/zm. The time for 
one scan was less than 4 hours, and the damage effects 
were negligible. Data count rate was about 10Hz. Ab- 
sorption is only about 6%, therefore it is neglected in the 
following analysis. 

7m(e _1 ) data of the Py-SO molecular crystal measured 
by IXS are shown in Fig.[2^i with momentum transfer q 
along the a axis. The spectral features are determined by 
the internal structure of the optically excited singlet ex- 
citons and interband transitions. There are three distinct 



features at about 2.2eV, 4.6eV and 6.6eV in the measured 
energy window, labeled as I, II, and III respectively. Af- 
ter removing the quasi-elastic Rayleigh background, the 
resulting spectra are shown in Fig.[2J>c. A weak feature 
IF and rising spectral weight beyond 8 eV can also be 
observed in Fig.[5fc. 

Exciton, or excited state in general has been a ma- 
jor challenge for computational physics or chemistry be- 
cause of the difficulty in treating electron correlation for 
excited states. Hutchison and coworkers have made a 
systematic investigation for the excited states of 60 or- 
ganic conjugated molecules with six commonly applied 
excited-state computational methods. When compared 
with experiments, ZINDO/SCI combined with Austin- 
Model- 1-optimized geometry was found to be the best 
choice in predicting the low- lying excited states, even out- 
performing the most commonly applied first-principles 
time-dependent density functional theory. [l|| Therefore 
it is adopted in our calculations. 

The IXS spectra calculated along the a axis are shown 
in Fig.[5Ji for a molecular aggregate of six-unit-cell stack 
with shortest intermolecular separations, [20j and a width 
of 0.17eV was included in order to compare with the 
experiment. One can also clearly identify several main 
spectral features A, B and B', C and D. There are al- 
most one to one correspondence between the experimen- 
tal features I-III and theoretical features A-C, respec- 
tively. The energy centroid positions of features B and C 
match those of features II and III almost perfectly. The 
weak feature B' at about 5eV becomes pronounced from 
1.12A -1 and above. The corresponding feature II' be- 
comes visible at high momentum transfers when feature 
II is weak. Both feature C and D involve many exci- 
tations. Correspondingly, the spectra in the experiment 
exhibit very broad feature III followed with rising spec- 
tral weight. These qualitative and, to a great extent, 
quantitative agreements between the theoretical spectra 
and the experimental results indicate that ZINDO/SCI 
calculations characterize the Py-SO system quite well. 

The experimental feature I disperses from 2.2eV at 
q = 0.28A- 1 to 2.07eV at q = 0.7A" 1 . This is well 
reproduced in the theoretical spectra, where the peak 
position of A disperses by 0.12eV in the corresponding 
momentum range, except that the theoretical position 
is about 0.48eV higher. The small exciton dispersion 
reflects the weak intermolecular coupling. In the 

aggregate calculation, the strongest intermolecular cou- 
pling between various orbitals in neighboring molecules 
is estimated to be 55meV. Therefore, the excitons are 
still Frenkel excitons that are confined mostly in a single 
molecule in this case. As a result, single molecule excita- 
tions could be computed to study the local distribution 
of the electrons and holes. In fact, S(q,u>) from single 
molecule calculations has been shown to agree very well 
with the EELS measurements for conjugated oligomers 
and polymers. [2l[ The IXS spectra calculated for molec- 
ular excitations of single Py-SO molecule along the a axis 
are shown in Fig. [2^ , where energy levels of the excited 
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FIG. 2: (a) IXS data of Py-SO with momentum transfer q along the a axis. The spectral intensities are divided by q 2 . (b, 
c) Enlargements of spectra shown in (a) with elastic line removed. ZINDO/SCI simulated S(q,uj)/q 2 for Py-SO for q along 
the a axis based on (d) the aggregate and (e) a single molecule. The dashed lines in panels a,c, and e indicate the calculated 
HOMO-LUMO gap energy. 



states are indicated by straight lines. As expected, the 
spectra are very similar to the aggregate calculations. 
The energy gap between LUMO and HOMO orbitals, is 
calculated to be 5.62 eV in ZINDO/SCI, denoted by the 
dashed line in Fig. [21 Moreover, our calculation shows 
that the lowest energy feature A corresponds to a discrete 
exciton, unlike in long chain systems, [12, HH where the 
peaks in the spectra might be consisted of many excitons. 
On the other hand, feature B contains mostly two ma- 
jor excitons and feature C and D are made up of tens of 
excitations above the gap. In experiment, as individual 
excitations could not be distinguished in features II and 
III, their dispersion could not be resolved even if there is 
any. 

As the dynamic structure factor is directly determined 
by the ground state and excited state wavefunctions, it 
puts stronger constraints on theoretical exciton distri- 
butions than just exciton energy position information 
derived from conventional optical measurements. The 
momentum dependence of the integrated weight is com- 
pared with the single molecule calculation in Fig.[3^-c 
for the three main features respectively. After rescaling, 
the theory agrees with the experiments very well. How- 
ever, as shown in Fig. 2, the theoretical high energy fea- 
tures are more pronounced than the experiment, produc- 
ing higher integrated weights in Fig. 3b-c compared to 
feature A. We attribute this overestimation to the more 
delocalized nature of the high energy excitons, which are 
less bound than the lowest exciton. Therefore, our sin- 
gle molecule calculation would overestimate its on-site 
occupation, which in turn causes higher matrix element 
of IXS. Indeed, the current aggregate calculation does 
give lower integrated weight, but more detailed calcula- 
tion on larger aggregate is needed to fully address this 
issue. As the large intermolecular distance corresponds 
to a momentum smaller than the sampled range, the cal- 
culated local distribution of the high energy excitons in 
a single molecule could still be compared with the ex- 
perimental data. For further confirmation, IXS spectra 
were sampled at two momentum transfers perpendicu- 
lar to the a direction in the ac plane, which also shows 



good agreement with the theory. (Fig. 4) The fact that 
the theory based on a single molecule matches the ex- 
periment so well indicates that the exciton distribution 
within the molecule is well captured, as a result of the 
weak coupling between the molecules. Fig.[5Ji-h are exci- 
ton wavefunction presented in a way that the false color 
scale indicates the possibility P{x, y) for finding an elec- 
tron at atom site x and a hole at atom site y. The gray 
scale of the solid circles in the inset shows the possibil- 
ity P(x) = J2 y P{ x i u) °f finding the electron or hole at 
site x on Py-SO molecule. These plots give clear infor- 
mation on how the electron and hole of a certain exciton 
are distributed in the molecule. For feature A, exciton 
is mostly situated in the middle region. For feature B, 
on the other hand, both of the two main excitons are ex- 
tended over the entire molecule, which are much larger 
than the lowest energy exciton. (Fig. 3e-f) For compari- 
son, the two main inter-gap excitations of feature III are 
shown in Fig.3g-h, one being confined in the phenyl ring, 
the other being extended. 

The measured intensity distribution of S(q,Lu) corre- 
sponds to the distribution of the exciton center-of-mass 
in the system. For simple linear molecular systems such 
as oligomer and polymer, S(q, uj) along the chain direc- 
tion relates directly with the exciton distribution over 
the molecule. The situation is more complicated for a 
molecular system with complex internal structures such 
as Py-SO. The structure of S(q, ui) measured along a 
certain q direction reflects distribution of the exciton 
center-of-mass in this direction summed across the en- 
tire molecule. Because the momentum transfer does not 
correspond to the relative electron-hole motion, quantum 
chemical calculations is needed to retrieve the full quan- 
titative details of the exciton wavefunctions. [13] Never- 
theless, there are still some qualitative correspondence. 
For example, the significant occupation of exciton A on 
atomic sites #15,16,19,20,21,24 (Fig. 3d) generally makes 
its center-of-mass more delocalized than the others along 
a-direction (Fig.l), and thus the momentum distribution 
of S(q, lo) is narrower in Fig. 3a. 

Inelastic x-ray scattering is a weak probe. From the ex- 
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FIG. 3: Integrated intensity as a function of momentum transfer for experimental features and single molecule calculations: (a) 
A and I integrated over [1.5eV, 3eV], (b) B+B' and II+IF integrated over [4eV, 5.4eV], and (c) C and III integrated over [5.4eV, 
8eV]. Broken and dashed lines in (b) and (c) show the q dependence of two main excitations of features B and C respectively. 
Fhe false color plots for the calculated possibility P(x, y) of finding an electron at site x and a hole at site y for (d) the lowest 
energy exciton; (e,f) two main excitations of feature B; (g,h) two main excitations of feature C. The atom sites are numbered 
in the top right corner. The atoms plotted in gray scale on the Py-SO molecule structure in the inset represent the probability 
P(x) = ^2 y P(x, y) of finding the electron or hole at site x for Py-SO molecule. 



perimental point of view, this is the first time that IXS is 
demonstrated to be feasible for organic molecular crys- 
tal in a third generation synchrotron. The experimental 
results are very clean and therefore can be directly com- 
pared with theory. The dispersion of the exciton feature 
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FIG. 4: (a) Experimental data with background removed. Mo- 
mentum transfer q is perpendicular to a axis, (b) ZINDO/SCI 
simulated S(q,ui)/q 2 based on a single molecule with q per- 
pendicular to a axis in the ac plane 



gives a good measure of the strength of the intermolecu- 
lar coupling. Combined with suitable quantum chemical 
calculations, reliable and comprehensive properties of ex- 
citons can be obtained, which are crucial for understand- 
ing their optical properties, and for designing materials 
of desired optical properties based on exciton transfer or 
dissociation properties. For example, if excitons in this 
useful energy range (2.2eV here) are localized within a 
particular complex structure, molecular clusters with a 
similar structure might be attached to polymers without 
affecting their local exciton (optical) behavior. 
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